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Abstract. It is shown that the number fx of free subgroups of index 6 A in the 
1 modular group PSL<2(Z), when considered modulo a prime power p a with p > 5, is 

always (ultimately) periodic. In fact, the main result of the paper is more general: 
| /' | 1 it establishes an analogous result for a one-parameter family of lifts of the modular 

. group containing PSL 2 (Z) as a special case. All this is achieved by explicitly de- 

I " termining Pade approximants to solutions of a certain family of Riccati differential 

equations. Our main result complements previous work by Kauers and the authors 
([Electron. J. Combin. 18(2) (2012), Article P37] and ["A method for determining the 
. mod-3 k behaviour of recursive sequences", preprint]), where it is shown, among other 

things, that the free subgroup numbers of PSL2(Z) and its lifts display rather complex 
■ behaviour modulo powers of 2 and 3. 

^— > ■ 

a- 
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1. Introduction 

T— I 1 

Beginning with the work of Sylow [25], Frobenius [5], [Bj, and P. Hall [9], [10], the study 
t^j- '. of congruences for subgroup numbers and related invariants has played an important 

Q\ ; role in group theory. Over the last fifteen years, there has been active interest in 

divisibility properties of subgroup numbers of (finitely generated) infinite groups; these 
may, to some extent, be seen as analogues of classical results for finite groups. The 
psj ! systematic study of the subgroup arithmetic of infinite groups begins with [17] . which 

investigates the parity of subgroup numbers and free subgroup numbers in arbitrary 
Hecke groups Sj(q) = C 2 *C q with q > 30 The results of [UJ were subsequently 
generalised to larger classes of groups and arbitrary prime modules in [I], [13], [IS] . 
US], [20], and [21]. 

A first attempt at obtaining congruences modulo higher prime powers was made in 
which studies the subgroup numbers of the inhomogeneous modular group Sj(3) = 
PSL 2 (Z) modulo 8, and derives a congruence modulo 16 for the number of free subgroups 
in PSL^Z) of given finite index. The real breakthrough concerning this kind of problem 
however occurred in the recent papers [IT] and [H], which introduce a method leading 
to semi-automatic existence proofs as well as explicit computation of such congruences 
modulo 2-powers and 3-powers for a wide variety of combinatorial sequences. Among 
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other things, this novel approach leads to congruences modulo arbitrary 2-powers for 
the number f\(m) of free subgroups of index 6mA in lifts T m of PSL^Z) of the form 

T m = C 2m * C 3m = (x,y\ x 2m = y im = 1, x 2 = y 3 ), m > 1; 

see Theorems 19 and 20 in [11]. The method developed in [11] can in principle be 
adapted to yield congruences modulo other prime powers p a , and this has been carried 
out in detail in [14J for p = 3. Among many other results, we obtain there congruences 
for the function f\(m) modulo arbitrary 3-powers; cf. Section 16 in [14J. The con- 
gruences obtained in this way for f\(m) show a highly non-trivial behaviour of these 
sequences modulo 2-powers and 3-powers. For instance, for the sequence f\ = f\(l) of 
free subgroup numbers of the group PSL 2 (Z), one can show that: 

(1) fx = —1 (mod 3) if, and only if, the 3-adic expansion of A is an element of 
{0,2}*1; 

(2) f\ = l (mod 3) if, and only if, the 3-adic expansion of A is an element of 

{0,2}*100* U {0,2}*122*; 

(3) for all other A, we have f\ = (mod 3); 

cf. [HI Cor. 53]. Here, for a set Q, we denote by Q* the free monoid generated by Q. 

The general theory of finitely generated virtually free groups implies that the gen- 
erating function for the free subgroup numbers f\(m) satisfies a Riccati differential 
equation, see (14.11) . This theory might also lead one to suspect that the primes 2 and 
3 should be special in this context. The main result of the present paper confirms this 
conjecture in a strong and surprising form. Namely, we show the following. 

Theorem. Let m be a positive integer, p a prime number with p > 5, and let a be a 
positive integer. Then the generating function F m (z) = 1 + X^a>i f^{ m ) zX i when coef- 
ficients are reduced modulo p a , can be represented as a rational function. Equivalently, 
the sequence {fx(m))x>i is ultimately periodic modulo p a . 

This is Theorem [H] in Section HI The proof relies on finding an explicit approximate 
rational solution for a relevant class of Riccati differential equations; see Theorem [T] in 
the next section. The proof of Theorem [1] occupies Section [3j Section H] presents the 
(very short) derivation of the theorem above from Theorem [TJ while Section [5] estab- 
lishes more precise results for the function f\ = f\(l) (i.e., the free subgroup numbers 
of the modular group PSL^Z) itself). More specifically, Theorem [11] identifies the 
denominators of the rational functions which one obtains when the coefficients of the 
generating function 1 + ^a>i /a 2 ^ are reduced modulo a prime power p a with p > 5. 
Theorems [T2TrT7l then illustrate this general result for the primes p = 7, 11, 13. In 
particular, the exact minimal period is determined in these cases. Finally, the rather 
interesting and involved way which led us to conjecture formulae (12. 3p and (12.41) for nu- 
merator and denominator of the approximate rational solution R n (z) of the differential 
equation (12. ip (which is the subject of Theorem [T]) is explained in an appendix. 

Note. This paper is accompanied by a Mathematics file that allows an interested 
reader to compute the rational function which provides the generating function 1 + 
Ea>i f° r t ne number f\ of free subgroups of index 6A in PSL 2 (Z) when the 
numbers are reduced modulo p a . The result is given as partial fraction expansion 
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of the form predicted by Theorem [TTJ The file is available at the article's website 
|http : //www. mat .univie . ac . at/~kratt/artikel/p sl2zmo d.htmTj 

2. Approximate rational solutions for a Riccati differential equation 

We consider the Riccati differential equation 

(1 - Az)F(z) - Bz 2 F'(z) - CzF 2 (z) - 1 - Dz = 0, (2.1) 

where A, B, C, D are constants. The technical main result of our paper is the following. 

Theorem 1. For every positive integer n, there exist uniquely determined polynomials 
P n (z) = 1 + Yl=\Vn,k zk and Qn{z) = 1 + Y2=i <ln,kZ k , where p n>k and q n;k are homo- 
geneous polynomials in A,B,C,D of degree k over the integers, such that the rational 
function R n {z) = P n (z)/Q n (z) satisfies the differential equation 

(1 - Az)R n (z) - Bz 2 R' n (z) - CzR 2 n (z) -l-Dz 

2n+l n 

= -t^tt(A + C + D) T\(£AB + AC + CD + i 2 B 2 + 2£BC + C 2 ). (2.2) 

Qui*) fJl 

Moreover, the coefficients p n ^ and q n ^ are given explicitly by 

( ]^n+l J^n—k 



Pn,n-k — „ fE 



2C{%-k) 2k+1 

k 

A+2C+E 



■(A + ^B+^E) ] (2.3) 



and 



n r>n—k 



■l) n B 



</<<■<<-<'• ~ ~ z IE _ ^ 



A+2C+E 



IB 

j=0 

k 



- E ( k+ k j ) (I _j) + i)^, ( M W s -) t J . (s-*) 

where E 2 = A 2 — ACD, and the Pochhammer symbol (a) m is defined by (a) m = 
a (a + 1)(ck + 2) • • • (a + m — 1) for m > 0, and (a) = 1. 
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Remark 2. (1) In the sums in (12. 3p it might seem that there is a problem if k = j — 
0. The correct way to interpret (12.31) for k = is the following: one first leaves k 
undetermined; then there is no problem to calculate the expression A + ^-B — jr^-E 
when j = 0, it simply equals A — E. Subsequently, one can safely set k = 0. A similar 
remark applies to its "companion expression" in the second sum. So, explicitly, 

Pn,n = { ~ ir 2 + CE n+l ~ E) - ^ + E ) ( M ^)n +1 ) " ( 2 ^) 

(2) The rational function solutions R n {z) for n— 1,2,3 are displayed in the appendix 
as {OMQ. 



(3) Formulae for the first few coefficients p n ^ and q^, k = 1,2,3, are displayed in 
(1A.4|) - (1A.9|) in the appendix. Assuming Theorem [H these (originally experimentally 



found) formulae can be established as follows: according to the theorem, the coeffi- 
cients p nt k and q ni k are in particular polynomials in A of degree at most k. As such, 
they are uniquely determined by k + 1 special evaluations. These can be found from 
Formulae (Q and fl23D by setting A = -2C - E - 2Bs, s = 0, 1, . . . , k. Indeed, 
with these specialisations, the term ( A+ 2^ +£ ) n+1 vanishes for n > k, so that the first 
terms between parentheses in (I2.3P and (12 .4p vanish. For the same reason, the term 
( A+ Ib +E ) ■ in the second 

sums over j in ( \2.3\\ and (12. Ah vanishes for j > k + 1, so that 
only the terms for j — 0, 1, . . . , k remain. 

(4) It is easy to see that R n (z) solves (12.21) if, and only if, R n (z/B) solves (12.21) with 
B replaced by 1, A replaced by A/B, C replaced by C/B, and D replaced by D/B. 
Hence, it suffices to consider the special case where B = 1. We shall occasionally make 
use of this simplification of the problem. 

(5) It is in fact not necessary to specify the precise form of the right-hand side in (12. 2p 
once we know that it is of the form 0(z 2n+1 ), that is, a formal power series of order at 
least 2n + 1; see Remark [51 Consequently, a different way to view the rational function 
R n {z) is as the Pade approximant of order n to the solution of the Riccati differential 
equation (12. ip . See also the proof of Lemma HI 



3. Proof of Theorem [Q 

The proof of Theorem [1] is accomplished by dividing the assertion of the theorem 
into several parts, and then treating each of these partial assertions in a separate 
lemma. More precisely, in Lemma [3] below we show that the rational function R n (z) = 
P n (z) /Q n ( z ) given by (12.31) and (12. 4p does indeed satisfy the differential equation (I2.2p . 
while uniqueness of solution is the subject of Lemma HI The fact that the coefficients 
p nt k and q n> k in (I2.3p and (I2.4p are homogeneous polynomials in A, B, C, D of degree 
k is then established in Lemma El while integrality of their coefficients is proved in 
Lemma [71 Altogether, these lemmas provide a full proof of Theorem [TJ 

Lemma 3. The rational function R n {z) = P n {z)/Q n {z), with the coefficients of the 
polynomials P n (z) andQ n (z) given by (\2.3h and (12.41) . satisfies the differential equation 



A RICCATI DIFFERENTIAL EQUATION AND FREE SUBGROUP NUMBERS 5 

Proof. If we substitute R n (z) = P n (z)/Q n (z) into the differential equation and multiply 
both sides by Qn(z), then we see that, in order to establish the lemma, we have to 
demonstrate the equation 

(I - Az)P n (z)Q n (z) - Bz 2 {P' n {z)Q n {z) - P n (z)Q' n (z)) - CzP 2 (z) - (I + Dz)Q 2 n (z) 

n 

= -z 2n+1 {A + C + D) Y[{£AB + AC + CD + i 2 B 2 + 2£BC + C 2 ). (3.1) 

We start by verifying that the coefficients of z 2n+l on both sides of (13.11) agree. For 
the sake of better readability, let us write 

n + := B^ e±§±^) n+1 and n_ := B^ (^) n+1 . (3.2) 

Using this notation plus (12. 5p . the coefficient p nj „ is seen to be given by 

Pn,n = {(A - E)U + - (A + E)IL) , (3.3) 

while, by (j2.4j) . the coefficient q n>n is given by 

(-l) n 

gn,n = ^-^(n + -n_). (3.4) 

Consequently, the coefficient of z 2n+1 on the left-hand side of (13 .ip equals 

-A Pn , n q n , n - Cp 2 n>n - Dql n = -^±^{{A{A - E)(2C) - C(A - E) 2 - D{2C) 2 )li 2 + 

+ ( - 2A 2 (2C) + 2C(A 2 - E 2 ) + 2D(2C) 2 )U + U_ 
+ (A(A + E){2C) - C(A + E) 2 - ,D(2L7) 2 )n^ 

= -^n + n_, 

where we have used the relation E 2 = A 2 — 4CD in the last step. It is straightforward 
to see that this is identical with the coefficient of z 2n+1 on the right-hand side of (13.11) . 

We now divide both sides of (13. ip by z 2n+1 , then differentiate them with respect to 
z, and finally multiply both sides of the resulting equation by z 2n+2 . In this way, we 
obtain the equation 

P n (z)(2CnzP n (z) - 2Cz 2 P' n {z) - (2n + l)Q n (z) + zQ' n (z) 

- z 2 {A + 2Bn - B)Q' n {z) + 2AnzQ n {z) + Bz 3 Q'^z)) 
+ Q n {z) (zP' n {z) -z 2 {A- 2Bn + B)P' n {z) - Bz 3 P^(z) 

+ (2n + l)Q n {z) + 2DnzQ n {z) - 2Dz 2 Q' n {z) - 2zQ' n {z)) = 0. (3.5) 

Since we already know that the coefficients of z 2n+l on both sides of (13. ip agree, es- 
tablishing (I3.5P is completely equivalent to establishing (13.11) . We concentrate on the 
former task from now on. 

The crucial observation is that, if (I3.5P is to hold, then the first factor between 
parentheses in (13. 5p must be a (polynomial) multiple of Q n (z), while the second factor 
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between parentheses must be a multiple of P n (z). Indeed, our claim is that 

2CnzP n (z) - 2Cz 2 P' n {z) - (2n + l)Q n (z) + zQ' n {z) - z 2 {A + 2Bn - B)Q' n (z) 

+ 2AnzQ n (z) + Bz 3 Q'^(z) = -{2n + 1 - nAz + n 2 Bz)Q n (z) (3.6) 

and that 

zP' n {z) - z 2 (A - 2Bn + B)P' n (z) - Bz 3 P' l [(z) + (2n + l)Q n (z) + 2DnzQ n (z) 

- 2Dz 2 Q' n {z) - 2zQ' n {z) = (2n + 1 - nAz + n 2 Bz)P n (z). (3.7) 

Clearly, if we manage to establish (13.61) and (13. 7p . then (13.51) . and hence also ( 13.1 ft and 
( \2.2\j follow immediately. 
Let us write 

Vn.k = pX k ^+ -p^, fe n_, 
q n ,k = <?„ ifc n + - g~ fc n_, 

where Il + and Il_ are given by ( 13.21) as before, which implicitly defines the p^ fc 's, etc., 
and let 

n n 

k=0 k=0 

n n 

Qn( z ) = QM = 

k=0 k=0 

One might hope that (13. 6p and (13 .7p continue to hold if one replaces P n (z) by Pn{z) 
and Q n (z) by Q^{z) everywhere, and that they continue to hold if one replaces P n (z) 
by P~(z) and Q n (z) by Q^{ z ) everywhere. This turns out to be too optimistic, but 
not too far from the truth. Our new claim is that 

2CnzP+(z) - 2Cz 2 (P+)'(z) - (2n + l)Q+(z) + z(Q+)'(z) -z 2 (A + 2Bn - B)(Q+)'(z) 
+ 2AnzQ+(z) + Bz\Q+)"{z) = -(2n + 1 - nAz + n 2 Bz)Q+(z), (3.8) 

2CWP-(^) - 2Cz 2 (p-)'(^) - (2n + 1)Q" (*) + z(Q~)'(z) - z 2 (A + 2Pn - B)(Q-)'(z) 
+ 2AnzQ-(z) + Bz 3 (Q-)"(z) = -{2n + 1 - nAz + n 2 Pz)Q;(;z), (3.9) 

z(P+)'(;*) - z 2 (A - 2Pn + B)(P+)'(z) - Bz 3 (P+)"(z) + (2n + l)Qt(z) + 2DnzQ + n (z) 
- 2Dz 2 (Q + J(z) - 2z(Q+)'(z) = (2n + 1 - nA* + n 2 &)P n + W 

2rA (-1)"5- (^§^L +1 



+(2 - +1) UJ cd^C/" ' (3 - 10) 

z(P-)'(*) - ^ 2 (A - 2Pn + B){P-)\z) - P^ 3 (P-)"(z) + (2n + l)Q~(z) + 2DnzQ~(z) 
- 2Dz 2 {Q-)'{z) - 2z{Q-J{z) = (2n + 1 - nAz + n 2 Bz)P~{z) 

2n\ (^i±^) n+1 
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These four equations do indeed imply (13. 6p and ( 13. 7ft : to be precise, by multiplying 
both sides of (13. 8p by H+ and multiplying both sides of (13. 9p by and then taking 
the difference of the two equations, we obtain ( 13. 6p . An analogous argument shows that 
( I3.10P and 03.1 lj) together imply 03. Tp . Moreover, since 

P n( z )\ E ^-E = -P~(z) and Q%(z)\ E ^_ E = -Q~(z), 

Equation (13.91) results from (13. 8p by replacing E by —E (note that the parameter 
D = (A 2 — E 2 )/(4C) remains invariant under this replacement), and an analogous 
assertion holds for ( I3.10P and (13. lip . Thus, the proof of the lemma has been reduced 
to establishing <g£E§ and ( 13TT0|) . 

We start with (13. 8p . Both sides of this equation are polynomials in z and A (with 
coefficients depending on n, B, C, E). Due to the form of the coefficients p^ k and q^ k , 
every subexpression in 03.8[) may conveniently (if tediously) be expanded in the basis 

It can be readily verified that corresponding coefficients on both sides of 03.81) agree. 

We use the same approach for (I3.10p . Again, one readily verifies that corresponding 
coefficients on both sides of 03 . 1 0[) agree, except possibly for the coefficients of basis 
elements of the form z° ( A+2 !^~ E ) - ; J ' — 0, 1, ... In other words, we have shown that 

( I3.10p holds except possibly for an error in the coefficient of z°. There are only three 
terms in (13.101) which have non-zero coefficient of z°: the term (2n + l)Q^(z) on the 
left-hand side, the term (2n + l)P^(z) on the right-hand side, and the last term on the 
right-hand side. So, in order to completely establish 03. lOj) . we have to show that 

(2n + 1)< - (2n + 1)< = (2n + 1) ( ™) { - , 2B )n+ \ 

or, explicitly, that 



■ (20 + A + %B - ^E) = (2» + 1) ^ n j ctg:^/"' - (3-12) 
The Gosper algorithm (cf. [7], § 5.7], [231 § II.5]) finds that 

.(2C + A + ^-B- gLE) = G(n,j + 1) - G(n,j), (3.13) 



with 



Clearly, summing both sides of 03 . 13[) over j running from to n immediately yields 
(1332]) . and thus OTTO]) . 

This completes the proof of the lemma. □ 



8 



C. KRATTENTHALER, AND T. W. MULLER 



Lemma 4. The differential equation (12. 2p has a unique solution R n (z) = P n (z)/Q n (z), 
where P n (z) and Q n (z) are polynomials in z of degree at most n, and Qn(0) = 1. 

Proof. The first observation is that the differential equation (12.11) has a unique formal 
power series solution F(z). This is seen by comparing coefficients of z m , m = 0, 1, . . ., 
which allows one to compute the coefficients of the series F(z) recursively. More pre- 
cisely, comparison of constant coefficients yields that -F(O) = 1, that is, that the coef- 
ficient of z° in F(z) equals 1, while, for m > 1, comparison of coefficients of z m yields 
an equation which expresses the coefficient of z m in F(z) uniquely in terms of lower 
(already computed) coefficients. 

Next one observes that the differential equations (12. ip and (I2.2p agree for the coef- 
ficients of z m , m = 0,1,..., 2n. Thus, what we are computing in (12. 2 p is the Pade 
approximant of the uniquely determined solution F(z) of (12. ip with the degrees of the 
numerator and denominator polynomials bounded above by n. Since Pade approxi- 
mants are unique (see e.g. [2]), a rational function R n (z) solving (12.21) . if there is one 
at all, is uniquely determined. However, from Lemma [3j we know that such a solution 
does indeed exist; and, by the arguments above, it is the unique solution of (12.21) . □ 

Remark 5. It should be observed that the above proof of Lemma @] did not use the 
precise form of the right-hand side of (12 .2p . Rather it shows that, if there is a rational 
function solution of (12. 2p of the prescribed form, then this right-hand side is forced 
upon us. 

Lemma 6. The coefficients p n ^ andq n ^, as defined in (I2.3P and (I2.4p . are homogeneous 
polynomials in A, B, C, D of degree k. In particular, p n $ = q n $ = 1. 

Proof. We first consider q n ^. We shall prove that the expression in (12. 4p is a polynomial 
in A, B, C, E. Taking into account that the denominator can be rewritten as 



and that the term between parentheses in (12 .4p is of the form Pol(-E) — Pol(— E), where 
Po\(E) is some polynomial in E, we see that q n ^ is a rational function in E 2 . Hence, if 
we are able to show that q n ^ is a polynomial in A, B, C, E, then, in view of the relation 
E 2 = A 2 — ACD, it will be obvious that q n ^ is also a polynomial in A, B,C, D. 

In order to accomplish this, we shall show that the expression between parentheses in 
( 12.41) vanishes for E = Bs with s = —k, —k + 1, . . . , k. This proves that this expression 
is divisible by the denominator on the right-hand side of (12. 4p (rewritten here in (I3.14p ) 
and, after multiplying numerator and denominator of the expression in (12.41) by B 2k+l , 
that q n ^ is a homogeneous polynomial in A, B, C, E of degree k, as desired. 

Writing the first term between parentheses, with E = Bs, in standard hypergeometric 
notation, 



k 




(3.14) 



i=i 
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we get 



' A+2C+Bs 
2B 



)n+l 

3=0 



fc + A (n — j 



(s+j + 1) fc ,( 



A+2C-B.S 



k — j J ' ' v 



' A+2C+Bs 
2B 



)n+l ( 



A+2C--Bs N 
2B 



+ s\ (n — s) 



(ra - fc)! 



X 3^2 



-fc + S, fc + S + 1, ~B 2 • ~y 

—n + s, s + 1 



. (3.15) 



On the other hand, if we write the second term between parentheses with E = Bs in 
hypergeometric notation, then we get 



( 



A+2C-Bs 
2B 



)n+l 



3=0 



( 



fc+j 
fc 



A+2C-Bs 
2B 



A+2C+,Bs 
2B 



), 



(fc + s) 



S! 



3-^2 



1* _|_ 1 .A + £ _|_ £ 
n,, ru -r i, 2B -r B -r 2 _ 

— n, s + 1 



(3.16) 



In order to see that one expression can be transformed into the other, we apply the 
transformation formula (see [U Ex. 7, p. 98]) 



3-^2 



a, b, c 
d,e 1 



T{e)T{d+e 



b-c) 



r(e-< 

6 = — fc, c 



T(d + e 
fc + 1, d 



"3^2 



a,d — b,d - 
d,d + e — b 



(3.17) 



with a =2s + ^ + f, & = —k, c = + d — s + 1, e — > —n to (13.161) (here, T(x) denotes 
the classical gamma function). After some manipulation, one sees that the result is 
exactly fl3TT5|) . 

The arguments for are completely analogous, albeit more involved due to exis- 
tence of the extra linear terms in A, C, E in the sums. One additional detail here is 
that one also has to show that C divides the expression between parentheses in (12. 3p . 
In order to see this, one sets C = in this expression. In view of E 2 = A 2 — ACD, this 
specialisation implies that E = ±A Without loss of generality, we let E — A. Then 
( 2B~ E )j = f° r a U J > 0- Hence, the only term possibly surviving is the one for 
j = in the first sum between parentheses in (12. 3p . However, for j = 0, C = 0, and 



A — tA = 0, so that this term vanishes also. 



E = Awe have A + ^B - ^E 

K+3 K+3 

Now, that we already know that p nj o and g nj o are constants, that is, that they are 
independent of A,B,C,D (or, equivalently, of A, B,C, E), in order to see that in fact 
Pn,o = <?n,o = 1, we are free to arbitrarily specialise the variables A, B, C, E. Our choice 



is A = — 2C — n — 1, B = 1, and E = —n — 1 . Then we have 



' A+2C-E 



2B 



). 



+ 1 



0. 



(0)„+i 

and therefore, in both (I2.3P and (12.41) . the second term between parentheses vanishes. 
For the same reason, in the remaining sums over j only the terms for j = survive, so 
that we obtain 

Pn,o = ^, ( ~S + l - l)n + i (» + 2)„ (-2C) = 1 



2C f-2n 



and, similarly, 



Qn,0 



l)2n+l 



(-2n - l) 2 „+i 



-n-l) n+1 (n + 2) B = l. 
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This completes the proof of the lemma. □ 

Lemma 7. The coefficients p n ^ andq n ^, as defined in (I2.3P and (j2.4|) . are polynomials 
in A, B, C, D with integer coefficients. 

Proof. By Lemma El we already know that p n ^ and q n ^ are polynomials in A, B, C, D. 
The assertion to be shown here is the integrality of the coefficients of these polynomials. 
In view of Part (4) of Remark [2l since we also know by Lemma [6] that p n ^ and q n ^ are 
homogeneous polynomials in A, B, C, D, we may without loss of generality set B = 1, 
which we shall do for the rest of the argument. 

By inspecting the expressions in (I2.3P and (12. 4p . one sees that the only possible 
obstacle to integrality of coefficients might be powers of 2 appearing in denominators. 

The claim is that no denominators occur after one has replaced E 2 by A 2 — 4CD. 
The reader should keep in mind that both expressions between parentheses in (12. 3p 
and (12 .4p are of the form Pol(E) — Pol(— E), where Pol(E) is some polynomial in E 
with coefficients depending on n, k, A, C . (We have again suppressed the dependence 
on variables other than E for the sake of better readability. The reader should recall 
that B — 1.) Moreover, the polynomial Pol(E) may be expanded as a polynomial in E 
and (A + 2C + E)/2. So, what we should examine are differences of the form 

E J A + 2C + E \ b JA + 2C-E' ' 



'A + 2C + E\ b [A + 2C 



Ea o ± — -S • ( 3 - 18 ) 




We claim that, after the substitution E 2 — > A 2 — ACD, expressions of this form have 
integer coefficients. Due to the sign ± in the previous expression, we would have to 
discuss two different cases. We content ourselves with treating the case where the sign 
is positive, the other case being completely analogous. We expand 

( A + 2C + E \ b ( A + 2C-E \" _ „ v 2"K , . h 

and then substitute A 2 — ACD for E 2 . Expanding again, this leads to the multiple sum 

2 2. £ 1 \£ 2 \(2£ 3 )\ A 6 U; 

We let a = l\ + 21^ d = £3 —£4. With this re-indexing, the above sum is converted into 
2 i-6 V y (-l) d 2 b -«b\(d + Q\ b _ a _ d d 

a^o Z&a ^ ~ 2 ^ ! (b-a-2d)\ (2d + 2£ 4 )\ £ A \ d\ 

a+d<b 

The sum over £4 can be written in hypergeometric notation, which leads to the expres- 
sion 

2 1^ AC D n n h _ n _oM(9M 2Fl 



a! (6- a -2d)! (2d)!' 

a+d<b 



a a I 1 

d+i - 1 
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This 2 -F\ -series can be summed by means of Gaufi' summation formula (see [24, (1.7.6); 
Appendix (III.3)]) 

T(c) r(c - a - b) 



2F1 

After some simplification, one arrives at 



a, b 

1 -L 



T(c- a)T(c-b)' 
D ^Td{ a ){b-a-2d) 



a,d>0 
a+d<b 



(-\) d A a C b - a - d D d ( ( 



a,d>0 
a+d<b 



a + d\ (a + d-i\\ ( b 
a J V a / / \b — a — 2d 



The integrality of coefficients in expressions of the form (13. 18)) is now obvious. As we 
argued above, this implies integrality of coefficients in q n ^, and, since — in comparison 
to q n> k — the expression (12.31) for p n ^ contains an additional 2 in the denominator, it 
implies that all coefficients of p n ^ are integers or half- integers (i.e., odd integers divided 
by 2). To see that the latter coefficients can in fact not be half- integers, let us suppose 
for a contradiction that p n ^ does contain half-integer coefficients, where k is minimal. 
Because of Lemma [6], which says in particular that p nj o = q n ,o = 1, we must have k > 0. 
We then consider the coefficient of z k on the left-hand side of ( 13. ip . It contains the 
term p n ,k<ln,o — Pn,k, and otherwise only terms involving p n j's with j < k and q n /s with 
j < k. Since the coefficient of z k on the right-hand side equals zero, this means that 
p U} k is given in terms of expressions which we know by the induction hypothesis to be 
polynomials in A, B, C, D with integer coefficients. The coefficient p n ^ must therefore 
have the same property, contradicting our assumption that p n ^ contains half-integer 
coefficients. This completes the proof of the lemma. □ 

4. Free subgroup numbers for lifts of PSL 2 (Z) modulo prime powers 

For integers m, A > 1, let f\(m) denote the number of free subgroups of index 6mA in 
the lift T m of PSL 2 (Z) defined by 

T m = C 2m * C 3m = (x,y\ x 2m = y 3m = 1, x 2 = y 3 ); 

(sun 

in particular, fx '■= /a(1) is the number of free subgroups in the inhomogeneous modular 
group PSIj2(Z) of index 6A. Our original motivation, when embarking on this project, 
was to establish the following result, which now is an easy consequence of Theorem [TJ 

Theorem 8. Let m be a positive integer, p a prime number with p > 5, and let a be a 
positive integer. Then the generating function F m (z) = 1 + 2~^a>i f\{ m ) zX j when coef- 
ficients are reduced modulo p a , can be represented as a rational function. Equivalently, 
the sequence (f\(m))\>i is ultimately periodic modulo p a . 

Proof. It is well-known that F m (z) satisfies the differential equation^ 

(1 - (6m - 2)z)F m (z) - 6mz 2 F' m (z) - zF^(z) - 1 - (1 - 6m + 5m 2 )z = 0. (4.1) 



2 Cf. [11, Formula (8.1)] for this, or Formula (18) in [H Sec. 2.3] with F rn (z) = 1 + 2/?r m(3) (z) for 
the same result in a more general context. 
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By specialising A = 6m — 2, B = 6m, C = 1, and D = 1 — 6m + 5m 2 in Theo- 
rem CE} we infer that there exist polynomials Pn{z) and Qn(z) such that .Ri^z) = 
Pt\z)/Q ( n l \z) satisfies 

(1 - (6m - 2)z)R^\z) - 6mz 2 (R^)' (z) - z {R^f (z) - 1 - (1 - 6m + hm 2 )z 

r. m 2n+2 2n+l n 

= -^1—4 TT(6^+ l)(6£m + 5). (4.2) 

(Qi m) ) w ir 

If we now choose n large enough then, since p a is coprime to 6, the right-hand side will 
vanish modulo p a , which proves our claim. □ 

5. Some more precise results on free subgroup numbers for PSL 2 (Z) 

In this section we shall concentrate on the number fx of free subgroups of index 6A 
in PSL^Z). Recall (see (14. ip with m = 1) that the generating function F(z) = 1 + 
Sa>i z x satisfies the Riccati differential equation 

(1 - Az)F(z) - 6z 2 F'(z) - zF 2 (z) -1=0, (5.1) 

which is the special case of (12. ip in which A — 4, B — 6, C — 1, and D — 0. 

The following two lemmas are the key for deriving the main result of this section, 
Theorem HU which identifies the denominators of the rational functions which one 
obtains when the coefficients of the generating function 1 + J2x>i f^ zX are reduced 
modulo a prime power p a with p > 5. In the statement of the lemmas, and also later, 
we write 

f(z) = g(z) modulo m 
to mean that the coefficients of z l in the power series f(z) and g(z) agree modulo m 
for all i. 

Lemma 9. Let p be a prime with p = 1 (mod 6). For A = A, B = 6, C = 1, D = 

(and, hence, E = 4), and n = ^Ezll ( mo d p), we have 

Qn{z) = Q( p -i)/ 6 (z) modulo p, 

with Qn(z) as given in the assertion of Theorem [B Furthermore, the polynomial 
Q{p-i)/q{. z ) has degree (p — l)/6 in z. 

Proof. We specialise the result of Theorem [1] to A = E = 4, B = 6, C = 1, D = 0, to see 
that R n (z) = P n (z)/Q n (z) solves (15.11) . where Q n (z) = J2k=o Qn,k zk , with coefficients 
given by 

^-PE((iueCI j ')G=9 ( - |+j+i) -^ 

-a +1 g( <: t i )6-3(i +j+1 ^^)- 

By Theorem [Q we know that q n $ = 1, and, from (13 .4p it is easy to see that 
<?(p-i)/6,0-i)/6 ^ (mod p) for A = E = 4, B = 6, C = 1, D = 0. (As a matter of fact, 
with this choice of parameters, we have Il + ^ (mod p), whereas Il_ = (mod p).) So, 
indeed, the polynomial Q( p _i)/ 6 (z) has degree (p — l)/6 in z. What remains to show is 
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that all coefficients q ntTl -ki < k < n — — 1, are divisible by p for the specialisation 
of the parameters A, B, C, D, E considered here. 

In order to see this, we start by writing the sums over j in (15.21) in hypergeometric 
notation, thus obtaining 



Qn,n—k 



-1 n 6 



nan—k 



2k+l 



V6/n+l 



|, k + 1, k 



3 ' 



—n 



;1 



(e)n 



n 

+ 1 u 



(I); 



3-^2 



g , fc + 1, /c 
5 ' 1 



3 - 



— n 



Next we apply the transformation formula (13.171) one more time, here with a = —k. 
This converts the last expression into 



•1P6 



nan—k 



Qn,n—k 



fc) 



2k+l 



V6/n+l 



■l) k 



+j A , (§)*(-»-*-§)* s*i> 

(i) (-n-k-%) , 3 F 2 

V6/n+l i. V3/feV 6/fc J z 



6' ^ 3' • 1 

i ir-t U 5 

3' ft 6 . 



2 —k 4- - —k 
6' h ^ 3' ■ ]_ 

5 ~i V. _1 

.3' " n 6 . 



;-i) n 6 



nan—k 



K (-1)*+* m (|-j) w+fc+ i _^ A /W\ (I ./')„ 

J/ U J/k+1 j=0 



E 

j=0 



fc! 



fc! 



(-i-i) 



fc+1 



(5.3) 



We claim that, if n = ^gi, 5 ^ ^ (mod p) and n > k + + 1, then each summand in 



the sums over j is divisible by p. 

To see this, we first consider the expression 

(|-j)n 



fc! 



•7' (3 ^)fc+l 



Let v p (a) denote the p-adic valuation of the integer a, that is, the maximal exponent e 
such that p e divides a. By variations of the well-known formula of Legendre [T6| p. 10] 
for the p-adic valuation of factorials, we have 



fc! 



^+fc+l 



(3 ^'/jfe+l 

p e 



E 



p £ 
1 ^ 



- j 



+ 



n + — j + l(p e + 5) 



fc-j + |G/ + 2) 



+ 



-J + |(^-1) 
7/ 



• (5-4) 



We need to prove that this sum over £ is at least 1. In the case where £ — 1, the 
summand reduces to 



i 




k-j 


.p. 




p 
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Let us write N = {n/p}, K = {k/p}, and J = {j/p}, where {a} := a — [a\ denotes 
the fractional part of a. With this notation, the last expression becomes 





n 




k 




.P. 




.P. 


+ 


N 







[K-J\ 

5 1 

h - 

6p 6 



6p ^ 6 



If n = (mod p), then N = Hence, 

N + K-J + — + - 
6p 6 



3p 3 



3p 3 



+ 



3p 



(5.5) 



Moreover, since k < n — — 1, we have 



> 1. Since, trivially, - [K - J\ > 



and 



~J L + 6 



+ 



J 3p + 3 



> 0, the expression in (j5.5p is at least 1. 



On the other hand, if n = 5 ^ 6 x ' (mod p), then N = 6v 

> 0, and, again, — —J 



. In this case, we have 



5 1 

N + K-J+ — + - 
bp o 



we have certainly 
If we suppose that 



n 




k 


V _ 




V _ 



6p ~ 6 



-J 



J_ 

3p 



> 0. 



•1, then K must be at least ^-{p — 1), or 



equivalently, k = 1), |(p— 1) + 1, . . . ,p— 1 (mod p). But, because of n > fc + £ g- + l 

> 1. Hence, in all cases the expression 



and n = 5 ^ ^ (mod p), this implies 



in ( 15. 5 p is at least 1. 

Now we have to discuss the summand in (15. 4p when £ > 2. Similarly to above, we 
write N = {n/p e }, K = {k/p e }, and J = {j/p e }- Using this notation, we may rewrite 
this summand as 



n 




k 


p l 







[K-J\ + 



N+K-J+-+ 



6 6 • p l 

K-J+- + 

3 3 • p l 



-J+-- 



+ 



6 6 • p l 
1 1_ 

+ 3 3-p*. 



. (5.6) 



We have - 
k + ^ + l, also 



J + 6 6/ 



+ 



J + 3 3y 



> 0, — [if — J\ ~> 0, and, since n > 



n 













> 0. If 



N+K-J+-+ 



6 6 • p e 



3 3 • p l 



> 0, 



(5.7) 



then (15. 6p is non-negative. There are two cases where (15.71) is violated. It should be 
noted that then the left-hand side of (I5.7P equals —1. Violation of (15. 7p can either 
happen if 

N + K - J + ~ + < < K - J + ~ + -^7 , 

O • p t 6 6 ■ p l 
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in which case necessarily K — J < and thus actually — \K — J\ = 1, or if 



N + K - J + - + ° 



1 2 

-. < 1 < K - J + - + j 

6 6 • p l 3 3 • pr 



In the latter case, we have in particular K > | — and on the other hand 



N < N + K — J 



3Y 

2 1 

< 



> 1, so that in all 



3 3 • p e 6 6 • p* 

Again, since n > k + + 1, this implies that actually 

cases the expression in (15.6j) is non-negative. 

This proves that the sum on the right-hand side of (15. 4R is at least 1, which, in turn, 
means that all summands in the first sum over j on the right-hand side of (15. 3ft are 
divisible by p. 

Next, we consider the expression 



k\ 



-k+l 



3 J H-3% +1 ' 



Here, we see that 
k\ 

oo 



( 3 ^)fe+l 

k-j 



J_ 

p ( 



pi 



+ 



rc + k- j + |(5 • j/ + 1) 



P 



p 1 



,j + l(2.j/-5) 



(5.8) 



Again, we have to prove that this sum over £ is at least 1. In the case where £ = 1, the 
summand reduces to 







k-j 


_p_ 




P 



+ 



n + k- j + i 



P 



+ 



6 



P 



+ 



6 



P 



+ 



2 




-i - 1 + 2 




+ 




3_ 


p 3 



With the notation N = {n/p}, K = {k/p}, and J = {j/p}, as before, the last expression 
becomes 



n 




k 


.P. 




.P. 



+ 



[K-J\ 



N + K-J+ — + - 
6p 6 



T 5 5 
-J h - 

6p 6 



2 2 
K-J + - 

3p 3 



+ 



7 5 2 

-J h - 

3p 3 



(5.9) 



If n = (mod p), then N = Hence, 

iV + A- J+ — + - 
6p 6 



La:- jj = i. 
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Moreover, since k < n — p -^- — 1, we have 
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> 1. If 



J- — + - 

6p 6 



+ 



n 




k 


P 




_P_ 






5 




J - 








3p 



then J - 2{p - 1] 2(p ~ 1} 
xnen j 3p , 3p 



,p — 1, and in all cases it follows that 

> 0. 



2 2 

K — J + - 

3p 3 



Therefore the expression in (15.91) is at least 1. 

On the other hand, if n = 5 ^~ 1 ^ (mod p), then A" = 5 ^~ 1 ^ . In this case, we have 



N + K - J+ — + - 
6p 6 



2 2 

K — J + - 

3p 3 



and 



3p 



n 








2(p-l) 



> 0. Again, if 



-J 



_5_ i 5 

6p 6 



J 



3p 



-1, then J 



3p 



+ 1, 



, p — 1 , and in all cases it follows that 

2(p-l 



- jj > 



3p 



> 0. 



If [If — J\ = (that is, if — [if — J\ ^ 1), so that in particular K > 



2(P~1) 
3p - 



then, since 



n 



5(P-1) 



(mod p) and n > + p -^- + 1, we actually have 



n 




k 


_P_ 




_P_ 



> 1. Hence, in all 



cases the expression in (15. 9p is at least 1. 

We turn to the summand in (15.81) when £ > 2. Writing A" = {n/p e }, K = {k/p e }, 
and J = {j/p e }, this summand acquires the form 



n 




k 









-J 



[K~J\ + 
5 5 



N + K -J- 
K-J 



6 6 • p 

Clearly, we have — \ K — J\ > 0, we have 

5 



N + K- J 



and, since n> k 



2=^ + 1, also 



1 

1 




"6 


• p e 




n 
pi 




k 

pi 




5 


+ 


6 


p l 



6 6 • p e 

2 2 

3 ~ 3-p e 



K-J 
> 0. If 



+ 



-J 



3-p e 



(5.10) 



3-p* 



> 0, 



7 2 5 

-JH 7 

3 3 V 



> 0, 



(5.11) 



then (I5.10P is non-negative. The inequality (15.111) is only violated if ~ — > J > 



6-p e 

gfj, in which case the left-hand side of (15. lip equals —1. It follows that K — 



J + 1 - < K + 4 < 1, so that 



3p e 
3 



K- J+Z-* 

3 3-p" 



< 0, that is, equivalently, 



> 0. If we now suppose that both — [K — J J and \N + K — J + 



A RICCATI DIFFERENTIAL EQUATION AND FREE SUBGROUP NUMBERS 



17 



I ^ 6^7 J are zero (that is, not 1), then K > J and N + K — J+| + < 1, and this 
entails the inequality chain 



l >N + K -J + l + -?—>N + K + - 
o • p p l 



On the other hand, since K > J > | — g^j, the number N can be at most 1/3. Because 

of n > + + 1, this implies that actually — -7 > 1, so that in all cases the 

expression in (I5.1(jp is non- negative. 

This proves that the sum on the right-hand side of (15. 8p is at least 1, which, in its 
turn, means that all summands in the second sum over j on the right-hand side of (15. 3ft 
are divisible by p, completing the proof of the lemma. □ 



Lemma 10. Let p be a prime with p = 5 (mod 6). For A = A, B = Q, C = 1, D = 

(and, hence, E = 4), and n = 2z_ ) -2^1 (mod p), we have 



Q n {z) = Q(p- B ye{z) modulo p, 



with Q n (z) as given in Theorem^ Furthermore, the polynomial Q ( p _5) iq(z) has degree 
(p — 5)/6 in z. 

Proof. This assertion may be established in a fashion similar to the proof of Lemma O 
The assertion about the degree of Q{ p -5)/q(z) is straightforward to see. 

In order to prove the congruence, we would again use the form (15. 3p for g„. n _fc. Here, 
we have 



(-i) fe+ V^ (I-j) 



' n+k+1 



£ 

1=1 

l even 



P £ 



k-j 



+ k-j + l(p e + 2) 



p' 



pi 



P 



+ 



+ 



n + k- j + l(p e + 5) 



-j + W-K 

p 1 



+ £ 



^ odd 





]_ 




k-j 


+ 


- 








pl 




pi 





n + k-j + \{f + 1) 

pi 



-j + \{hp l -l) 



pi 



k-j + l(p e +l) 
pi 



+ 



-j + 

pi 



(5.12) 
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and 



i) k+j fk\ a-j) n 



+k+l 



oo 



1=1 

I even 

j + |(p f -l) 

pi 









k-j 


+ 


- 


p^ 




p e 



n + fc- j + i(5p* + l) 



fc-j + |(p f -l) 



+ 



-j + §(2p*-5) 



£ odd 







k-j 


+ 










pi 





n + k-j + + 1) 

pi 



-J + |(^-5) 
pi 



k-j + ±(2^-1) 

pi 



+ 



-j + I(2/-4) 



(5.13) 



and these are the expressions which have to be analysed p-adically. We leave it to the 
interested reader to fill in the details. □ 

Theorem 11. Let a be a positive integer. If p is a prime with p = 1 (mod 6), then 
the generating function F(z) = 1 + X1a>i/a-2 A ? when coefficients are reduced mod- 



ulo p a , equals P a (z)/Q 



(p-l)/6 



[z), where P Q (z) is a polynomial in z over the integers, 
and Q( p -\)/q(z) is the polynomial of degree (p — l)/6 given by the formula in Theorem^] 
with A = E = 4, B = 6, C = 1, and D = 0. 

On the other hand, if p is a prime with p = 5 (mod 6), then the generating 
function F(z) = 1 + Y2\>i f^ zX > when coefficients are reduced modulo p a , equals 
Pa{ z )/Q( p -5)/e( z )7 where P a (z) is a polynomial in z over the integers, and Q( p -5)/6(z) is 
the polynomial of degree (p — 5)/6 given by the formula in Theorem^ with A = E = 4, 
B = 6, C = 1, andD = 0. 

Proof. We proceed by induction on a. 

For the start of the induction, we just have to use Theorem [T] with A = E = 4, 
B = 6, C = 1, D = 0, and choose n = (p — l)/6 and n — (p — 5)/6, respectively, as 
this makes the product on the right-hand side of (I2.2p vanish. By uniqueness of the 
solution of (15.11) . it follows that 



F(z) 



Pd(z) 
Q d {z) 



modulo p, 



where d = (p— l)/6 or d = (p — 5)/6, depending on the congruence class of p modulo 6. 

Let us now suppose that we have already found an integer polynomial A a (z) such 
that F a (z) := A a (z)/Q2(z) solves (15. ip modulo p a for some a > 1. We choose n = 
(mod p) or n = 2^ (mod p), respectively, large enough such that the product over 
£ on the right-hand side of (14 .2p with m = 1 vanishes modulo p a+l . (The reader should 
recall that the specialisation A — E — 4, B — 6, C — 1, D = 0, which is considered here, 
corresponds to the specialisation of Theorem [T] considered in Theorem [H] for m = 1). 
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By uniqueness of solution of the differential equation (14. 2p (see Theorem [TJ , we must 
then have 

R n (z) = F a (z) modulo p a , (5.14) 

where R n (z) = R [ n\z) = P n (z)/Q n (z) is the solution to (14.21) with m — 1 given by 
Theorem [TJ Consequently, if we consider the difference, 

p M _ p M - IjM. _ i«(f) - ^n(*)Qg(*)-4,(s)Q n (s) 

then, since Qg(2)Q Tl (z) has constant coefficient 1 (see Theorem [TJ), we know by (I5.14p 
that all coefficients of the integer polynomial in the numerator of the last fraction must 
be divisible by p a . In other words, there is an integer polynomial B a (z) such that 

B a (z) 



a 



If we consider both sides of this equation modulo p a+l , then the fraction which is 
multiplied by p a may be reduced modulo p. By Lemmas M and [TUl this leads to the 
congruence 



A a (z) a B a (z 



Mz) = : ^^+P a Zl\~\ modulo p a+1 . (5.15) 



z 



So, indeed, the solution R n (z) to the Riccati differential equation (15.11) . when coeffi- 
cients of series are reduced modulo p a+1 , can be expressed as a rational function with 
denominator Q^ +1 (z), where d — (j> — l)/6 in case p = 1 (mod 6), and d = (p — 5)/6 if 
p = 5 (mod 6). This establishes the theorem. □ 

In what follows, we shall make Theorem [11] more explicit for the first few prime 
numbers, and determine precisely the (minimal) period length. 

Theorem 12. Let a be a positive integer. Then the generating function F(z) = 1 + 
fx z x , when coefficients are reduced modulo 7 a , equals P a (z)/(1 + 2z) a , where 
P a {z) is a polynomial in z over the integers. In particular, the sequence (/a)a>i is 
ultimately periodic, with minimal period equal to 6 • 7 a ~ l . 

Proof. In view of Theorem [TTl and since Qi(z) = l + 2z modulo 7 for A = E = 4, B = 6, 
C — 1, the only assertion which still requires proof is the one about the period length. 
A careful examination of the previous arguments (see in particular (15.151) ) reveals that 
we have actually shown the stronger statement that the solution F(z) to (15. ip . when 
coefficients are reduced modulo 7 a , can be written in the form 

a 

F(z) = Po\ a (z) + + 2z)~ fc modulo 7°, 

k=l 

where Pol a (z) is a polynomial in z with integer coefficients, and the /_jt's are integers. 
Now, we have 



7 fe - 1 (l + 2z)- k = ^(-2) £ 7 fc - 1 ( k + k l _ x l \ 



£>0 



The term (— 2) , £ = 0, 1, . . ., is periodic modulo 7 a with (minimal) period length 
(p(7 a ) = 6 ■ 7 Q_1 , where (p denotes Euler's totient function. Indeed, this holds for 
a = 1, 2 by inspection, and thus for every a > 1 by a theorem of Gaufi concerning 
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existence of primitive roots (cf. e.g. [3^ Ch. 2, Sec. 5] or [221 pp. 285-287]). On the 
other hand, the period length of 

7 k-i (k+e-i\ k ^_ M{k _ m {k + ...(i + 2)^ + 1) 

V fc-1 J (k - i)\7-M(k-m ' " ' 

when considered modulo 7 a , divides 7 a_1 for k > 2, since 

k- 1 -v 7 ({k- 1)!) > 1 

for k > 2. In total, 

(- 2 )' 7fe_1 ( /c t-T 1 )' f=0 ' 1 '---' 

when considered modulo 7 a , is periodic with period length (exactly) 6 • 7 a ~ 1 , which 
implies the corresponding assertion of the theorem. □ 

We have implemented the algorithm which is implicit in the proof of Theorem [Til (see 
the Note at the end of the Introduction). As an illustration, the next theorem contains 
the result for the modulus 7 5 = 168 07. 

Theorem 13. We have 

A>1 

= 4802z 25 + 9604z 23 + 14406z 22 + 2401z 21 + 2401z 20 + 4802z 19 + 9947z 18 + 9604z 17 
+ 1029(k 16 + 99472 15 + 10976^ 14 + 16464z 13 + 12691z 12 + 2940Z 11 + 8918z 10 
+ 15484z 9 + 8722z 8 + 4214z 7 + 10829/ + 6174/ + 406z 4 + 14896z 3 + 11102z 2 

„ 16451 9562 2450 2744 2401 
+ 14168^ + 7+ — + — r + — - + - — - + 



l + 2z (l + 2z) 2 (l + 2z) 3 (l + 2z) i (l + 2^) 5 

modulo 7 5 . (5.16) 

The case of powers of 11 can be handled in an analogous way. We content ourselves 
with stating the corresponding result, leaving the proof to the interested reader. 

Theorem 14. Let a be a positive integer. Then the generating function F(z) = 1 + 
J2x>i /a z X , when coefficients are reduced modulo ll a , equals P a (z)/(1 — z) a , where 
P a (z) is a polynomial in z over the integers. Moreover, the sequence (fx)x>i is ultimately 
periodic, with minimal period length equal to ll a_1 . 

Again, we have implemented the algorithm which leads to the above theorem (see 
the Note at the end of the Introduction). As an illustration, the next theorem contains 
the result for the modulus ll 5 = 161051. 

Theorem 15. We have 

A>1 

= 87846/ 1 + 87846/ 9 + 131769/ 8 + 87846z 37 + 146410/ 6 + 29282/ 5 + 87846/ 4 
+ 87846^ 33 + 131769/ 2 + 123783z 30 + 146410/ 9 + 65219/ 8 + 151734/ 7 
+ 153065/ 6 + 105149z 25 + 154396/ 4 + 145079/ 3 + 153065/ 2 + 22627/ 1 
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+ 1038182 20 + 4719/ 9 + 78529/ 8 + 156453/ 7 + 153186/ 6 + 64614/ 5 
+ 1231782 14 + 20933z 13 + 154033/ 2 + 84579/ 1 + 93533/° + 1514922 9 
+ 283252 s + 136730/ + 23727/ + 43164/ + 75636/ + 1493582 3 + 126445/ 

„ 80547 6809 17787 41261 14641 
+ 973832 + 7 + + - + - + - + 



, . (I-2) 2 (I-2) 3 (l-^) 4 (l-^) 5 

modulo ll 5 . (5.17) 

The next prime, 13, is congruent to 5 modulo 6, and Theorem [TT] says that we 
may present solutions modulo 13 e as rational function with denominator being the e-th 
power of a quadratic factor. It turns out that this quadratic factor is (1 — 2z)(l + 5z). 
Experimentally, it seems that the maximal exponent s such that (1 + 5z)~ s appears 
in the (cleared) denominator is actually s — (e + l)/2 (instead of e; see for example 
Theorem [TTj) . but our methods are not able to demonstrate this. What we are able to 
prove is summarised in the theorem below. 

Theorem 16. Let a be a positive integer. Then the generating function F(z) = 1 + 
^\>i f\z x , when coefficients are reduced modulo 13 a ; equals P a (z)/ ((1 — 2z)(l + 5z)) , 
where P a {z) is a polynomial in z over the integers. Moreover, the sequence {fx)x>i ^ 
ultimately periodic, with minimal period length equal to 12 • 13 a_1 . 

The algorithm which led to the above theorem generates the following result for the 
modulus 13 5 = 371293. 

Theorem 17. We have 

A>1 

= 314171z 42 + 2856102 40 + 1428052 38 + 114244z 37 + 2856102 36 + 571222 35 
+ 118638z 34 + 2856102 33 + 325156z 32 + 142805z 30 + 90077z 29 + 3383382 28 
+ 3493232 27 + 1889422 26 + 103259z 25 + 26364z 24 + 35 1 52z 23 + 1889422 22 

+ 4732z 21 + 76895z 20 + 3106222 19 + 285612 18 + 340535z 17 + 3587872 16 
+ 3533792 15 + 1350312 14 + 1155962 13 + 20280z 12 + 3288742 11 + 559392 10 
+ 1164412 9 + 56745/ + 179309/ + 342212/ + 219700/ + 24336/ 
+ 238953/ + 332462/ + 354965;? + 13 
208033 363181 171366 334822 
+ l + 5z + (1 + 5z) 2 + (1 + hzf + 1 - 2z 
176228 154635 134017 314171 , , - , 

+ 7 777 + 7 7T + 7 77 + 7 77T modulo 13 5 . (5.18) 

(1-22) 2 (1-22) 3 (1 - 22) 4 (1 - 22) 5 V ' 

For p = 17, the denominators of rational solutions to ( 15. ip (viewed modulo a power 
of p) are powers of 1 + 152 + 7/ (which does not factor modulo 17), and for p = 19, 23 
they are powers of (irreducible) ternary factors, etc. A generalisation of the argument 
in the proof of Theorem [12] shows that it is still true that the period of fx modulo a 
prime power p a is a multiple of but, which multiple it actually is, is in general 

difficult to describe. For example, modulo 17 the period of fx is 6 ■ 16, modulo 17 2 it is 
18 ■ 16 ■ 17, while modulo 17 3 it is 102 • 16 • 17 2 . 
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Appendix: How were the expressions for p U)k and g„ )fc found? 

When papers get written, the path(s) towards the results presented often get lost, are 
covered up, or — in some cases — even hidden intentionally. As we believe that, in the 
present case, this path is actually quite interesting and instructive, we shall describe in 
this appendix how we were led to conjecture the formulae in (12. 3p and (12.41) . without 
which we would not even have been able to prove the mere existence of solutions to 
(I2.2p . and without which we would therefore never have reached the result described 
in Theorem [H] concerning the number of free subgroups in lifts of the inhomogeneous 
modular group, which was our primary goal and motivation. 

Recall that the generating function F(z) = 1 + J2x>i z> "i wnere fx denotes again 
the number of free subgroups of PSL^Z) of index 6A, satisfies the Riccati differential 
equation (15. ip . This equation allows us to compute the numbers fx recursively. We 
started out by looking at tables of these numbers, when they are reduced modulo powers 
of 5, 7, 11, 13, etc. The case of 5 is a trivial (and exceptional) one, since it is not hard 
to see that f\ vanishes modulo any power of 5 once A is large enough. The next case, 
powers of 7, is more interesting. Computer experiments quickly led to the conjecture 
that the free subgroup numbers fx modulo a given power of 7 are ultimately periodic, 
and it was even possible to predict the period length. Theorem [12] presents the precise 
corresponding statement. Continuing with powers of 11, a similar picture emerged. 
The precise statement is presented in Theorem [141 With some effort, we could see that 
the free subgroup numbers f\ are also periodic when reduced modulo 13 or 17, however 
with considerably larger period lengths. What happens modulo higher powers of 13 or 
17 was unclear at that point. 

We then decided to try to prove, say, the statements modulo powers of 7 and 11. An 
approach in the style of [11] and [H] came quite close, but did not succeed. At some 
point, we had the idea of "Pade- approximating" the solution to (15.11) by a rational func- 
tion with numerator and denominator degrees equal to n, and see how far this rational 
function fails to satisfy the differential equation ( 15. ip . We were (pleasantly) surprised 
to see that the error on the right-hand side was a "round" | multiple of z 2n+1 /Q 2 n (z), 
and, using the guessing programme Rate3 it was not difficult to come up with the 
conjecture that (14. 2 p (with m = 1) should hold. 

After several failed attempts at proving this conjecture, in a mood of despair, we de- 
cided to make everything more general, not expecting anything but horrendous results. 
We put as many parameters as possible into the differential equation (15.11) . so that we 
were led to consider the differential equation (12.11) . We were in for another pleasant 
surprise: although one cannot go very far with nowadays' computers (we only made it 
up to n = 5), the pattern which is displayed in what is now Theorem [1] became quickly 
clear, however without any clue about the explicit form of the rational function R n (z). 



3 According to Kuperberg (T51 Sec. VILA], "an integer is round if it is a product of relatively 
small (prime) numbers. A round enumeration in combinatorics almost always has an explicit product 
formula." See [TH Sec. 2] for a more detailed elaboration of the "notion" of "roundness." 
4 written in Mathematica by the first author; available from 
: //www . mat . univie . ac . at/~kratt 
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For example, we obtained 

r> , \ l + (-B-C + D)z ^ 

R M = l + {-A-B-2C)z > (A ' 1} 

1 + {-A - AB - 3C + D)z + (—AD + 2B 2 + 3BC - WD + C 2 - 3CD) z 2 

2 ^ Z > ~ 1 + (-2A -AB- AC)z + (A 2 + 3AB + 3 AC + 2fi 2 + 6BC + 3C 2 - CD) z 2 ' 

(A.2) 

* w = (A - 3) 

with 

P 3 (*) = 1 + s(-24 - 9B - 5C7 + £>) 

+ z 2 (A 2 + 7AB + AAC - 2AD + 18B 2 + - 8RD + QC 2 - QCD) 

+ z 3 (A 2 D + 6ABD + AACD - QB 3 - \\B 2 C + \\B 2 D 

- 6BC 2 + 18BCD -C 3 + 6C 2 D - CD 2 ) 

and 

Q 3 (z) = 1 + z(-3A -9B- 6C) 

+ ^ 2 (3A 2 + 15AB + KMC + 185 2 + + IOC 2 - 2CD) 

+ z 3 (-A 3 - 6A 2 B - AA 2 C - MAB 2 - 18ABC - 6AC 2 + 2ACD - 6B 3 

- 22B 2 C - 18BC 2 + 6BCD - AC 3 + AC 2 D). 

On the basis of these computer data (including the ones for R±{z) and R 5 (z)), it is not 
difficult to come up with guesses for the first few coefficients p Ht k and q n ^: 

p n l = -{n - I) A - n 2 B - (2n - l)C + D, (A.4) 

q n l = - n A - n 2 B - 2nC, (A.5) 

Pn,2 = \{n - 2)(n - l)A 2 + \(n - 2)(ra - l)(2ra + 1)AB + 2(ra - 2)(n - 1)AC 

-{n- I) AD + |(n - l) 2 n 2 B 2 + (n - 1) (2n 2 - 2n - l) fiC 

- (n - l)(n + + {n - l)(2n - 3)C7 2 - 3(n - 1)CT>, (A.6) 
5n,2 = - 1M 2 + \(n - l)ra(2n - 1)AB + (n - l)(2n - 1)AC + |(n - 1 ) 2 n 2 B 2 

+ {n- l)n{2n - l)BC + (n - l)(2n - 1 )C7 2 - (n - 1)CD, (A.7) 
p n3 = -I( n - 3)(n - 2)(n - f)A 3 - ±(n - 3)(n - 2) (n 2 - ra - l) A 2 B 

- \{n - 3){n - 2){2n - 3)A 2 C + \{n - 2)(n - f)A 2 D 

- J(n - 3)(n - 2) (3n 3 - 3n 2 - n - 2) AB 2 

- \{n - 3){n - 2)(2n - 3){2n + l)ABC + ±(n - 2)(n + l)(2n - 3)ARD 

- (n - 3)(ra - 2)(2n - 3) AC 2 + (n - 2)(3n - h)ACD 

- \{n - 2) 2 {n - l) 2 n 2 B 3 - \{n - 2){2n - 3) (3n 3 - Qn 2 - n - 2) B 2 C 
+ \{n - 2) {n 3 - n - 2) B 2 D - (n - 2)(2n - 3) (n 2 - 2n - l) BC 2 

+ (n - 2) (3n 2 - 2n - 3) BCD - \{n - 2){2n - 5)(2n - 3)C 3 



24 



C. KRATTENTHALER, AND T. W. MULLER 



+ 2{n-2){2n-3)DC 2 - (n-2)CD 2 , (A.8) 
q n 3 = -l( n - 2)(n - l)nA 3 - \(n - 2)(n - l) 2 nA 2 B - (n - 2)(n - 1) 2 A 2 C 

- \{n - 2){n - l)n (3n 2 - 6n + 2) AB 2 - (n - 2){n - l)n(2n - 3)ABC 
-(n-2)(n- \){2n - 3) AC 2 + (n - 2)(n - 1)ACD 

- \{n - 2) 2 {n - l) 2 n 2 B 3 - |(n - 2)(n - l)n (3n 2 - 6n + 2) S 2 C 

- (n - 2)(n - l)n(2n - 3)5C 2 + (n - 2)(n - l)nBCD 

- |(n - 2)(n - l)(2n - 3)C 3 + 2(n - 2)(n - l)^ 2 ^. (A.9) 

So, it seemed "clear" that the p n ,kS and q n ,k s were homogeneous polynomials in A, B, C, 
D of degree k. However, having stared at this for some while, it seemed hopeless to us 
to continue along these lines. 

We made the "intermediate" observation that, should we be able to demonstrate 
something at all in this context, it would have to come from the "hypergeometric world," 
and consequently the quadratic factors in the right-hand side product of (12.21) must 
factor into linear factors. This happens only if the discriminant of the factor, viewed 
as a quadratic form in £, has a square root, and the latter requires the parametrisation 
E 2 = A 2 — 4CD. If one uses this parametrisation, then the right-hand side product in 
(12. 2p becomes 

n 

-{A + C + D) Y[(£AB + AC + CD + £ 2 B 2 + 2£BC + C 2 ) = --n+n_, 

i=i 

where, again, we make use of the short-hand notation in (13.21 ). 

At a certain point, we decided to "look at the other end," that is, to try to find 
formulae for the coefficients p n ^ and q n ^ for k = n, n — 1, . . . Comparing coefficients of 
z 2n+i on both sides f (|3. ip . one obtains the equation 

-Ap ntn q n , n - Cp 2 nn - Dq 2 nn = -^n + n„. 

Moreover, upon using the parametrisation E 2 = A 2 — 4CD, the left-hand side factors, 
leading to 

(Cp n , n + \(A- E)q n:n ) (Cp n , n + \{A + E)q n>n ) = n+IL. 

If there is any "justice" in the world, then one of the factors on the left-hand side would 
have to agree with one of the factors on the right-hand side, and the same should hold 
for the other two factors. Indeed, the computer told us that, apparently, 

Cp njn + \{A- E)q n>n = n_, 

Cp n:n + l(A + E)q nin = U+. 
Solving this system of linear equations for p n ^ n and g n>n , one arrives at 

Pn,n = K ~x^(( A - £) n + ~( A + E)n-) 



(that is, Equation ( 13. 3p ). and 



(that is, Equation ( l3.4p ) 



(-l) n 

q n ,n = ^M n + " n - 
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Subsequently, we proceeded to compare coefficients of z 2n on both sides of (13. ip . 
This leads to a linear equation for p n ^ n -\ and q n , n -\i with coefficients containing U + 
and n_. Containing two unknowns, it is certainly undetermined. Nevertheless, one is 
able to come up with a guess: one considers the equation first modulo Il + — and 
then modulo Il + + As it turns out, this leads to dramatic simplification, and one 
obtains two linear congruences for p n ,n-i and q n ,n-i (the moduli being Il + — n_ and 
n + + respectively), with coefficients being rational functions in A, B,C, E. One 
can then again use computer experiments to come up with a guess for the multiplicative 
factors which make these congruences into equations. Once this is done, one has two 
linear equations, which one solves for p n>n _i and q n ^ n -i- The result was ( 12. 3 p and (12 .4p 
with k = 1 (albeit not yet written in this form). We did the same for p n>n -2 and g ni „_2, 
and for p nj „_3 and q n , n -3- (The latter required to also compute R§(z) and Rj(z), which 
now, by putting the knowledge of the coefficients p n ,k and q n ^ for k = n, n — 1, n — 2 
into the computer programme, became feasible.) 

At this stage, we realised that the coefficients of H+ and Il_ became "elegant," if one 
expanded them — viewed as polynomials in n — in the basis 

1, (n-k + 1), (n - k + \){n - k + 2), (n - k + l)(n - k + 2) • • • (n - l)n. 

We also had at this stage abundant "information" on the arising denominators and 
signs. More precisely, our experimental data told us that one could apparently write 

^ l) n + 1 J^n—k 

Pn,n—k 



2C(f -k) 2t+1 

j=0 

-{ A±2 ^) n+ ^- k + l )^(-B+^ l )k- ] \ 



k 

A+2C+E \ 
2B )j 
3=0 



{P { uIA + PnL B + PnL E ) 



and 



-l) n B 



n r>n—k 



qn,n-k ~ , E n 

Kb ~~ K )2k+1 

( i± W s )„» t **,(» - k + (-1 + i + 

3=0 

3=0 

(A.10) 

with "nice" numerical coefficients Pn\ j^Pni j^Pnk j^n,k,j- The available data at this 
point (namely, the rational functions R n (z) for n = 1,2, ... ,7) were largely insufficient 
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to make any guesses about these coefficients. However, assuming the above forms of p n ^ 
and <7 n ,fc, together with the earlier observed polynomiality of p n ,k and q n ^ in A, B, C, D, 
and, hence, also in A, B, C, E (except for possible powers of C in denominators), suffi- 
cient data could now be easily produced! Let us first concentrate on q n , n -k (as we did). 
If q n , n -k is to be a polynomial, then (-| — k) 2 k+i must divide the expression between 
parentheses on the right-hand side of (1A.10I) . that is, this expression must vanish for 
E = Bs, s = —k, —k + l,...,k. Let us do the substitution E = Bs in this expression, 
and let us suppose that s > 0. Then we have 

( A+2C+E \ ( A+2C+Bs \ _ ( A±2C±Bs\ ( A+2C-Bs , , i\ 

I 2B )n+l \ 2B ) n +l V 2B ) n+ i- s \ 2B ~ rU ^~ 1 ) s ' 

Comparing with 

I A+2C—E \ __ ( A+2C-Bs \ 

\ 2B )n+l \ 2B /n+1' 

we infer that ( A+2 £b~ Bs + n + l) must divide the second sum over j as a polynomial 
in n. This provides many non-trivial vanishing conditions for this sum over j, viewed 
as polynomial in n (they actually "overdetermine" it), and this suffices to compute the 
coefficients q n ,k,j f° r a large range of fc's and corresponding j's. Using the guessing 
programme Rate (see Footnote H|) once again, it is then straightforward to guess the 

formula in (12. 4p . For the coefficients Pn\ jiPn\ j^Pn\ j one d° es n °t have to do the same 
work again: it suffices to look at the quotients 

(1) (2) (3) 

"n,k,j "n,k,j "n,k,j 

5 5 5 

Qn,k,j Qn,k,j Qn,k,j 

which quickly leads one to the formula in (12.31) . 
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